A theoretical study of tunneling conductance in PrOs 4 Sbi 2 superconducting junctions 
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The tunnel conductance in normal-metal / insulator / PrOs4Sbi2 junctions is theoretically stud- 
ied, where skutterudite PrOs4Sbi2 is considered to be an unconventional superconductor. The 
conductance is calculated for several pair potentials which have been proposed in recent work. The 
results show that the conductance is sensitive to the relation between the direction of electric cur- 
rents and the position of point nodes. The conductance spectra often deviate from the shape of 
bulk density of states. The sub gap spectra have peak structures in the case of the spin-triplet pair 
potentials. The results indicate that the tunnel conductance is a useful tool to obtain information 
of the pairing symmetry. 

PACS numbers: 74.50.+r, 74.25.Fy,74.70.Tx 
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I. INTRODUCTION 

Superconductivity in the cubic skutterudite PrOs4Sbi2 
(POS) has received a great interest in recent years since 
it has two superconducting phases 1 . The specific heat 
results 2 show jumps at T c \ = 1.82K and T c2 = 1.75 
K. Nowadays such two superconducting phases are well 
known in a spin-triplet superconductor UPt3 and a su- 
pcrfluid 3 He. A NQR experiment shows the absence of 
the coherence peak, which suggests that POS is an un- 
conventional superconductor—. A thermal conductivity 
experiment indicates 6 point nodes at (1,0,0) direction 
and directions equivalent to (1,0,0) for the high temper- 
ature phase (H-phase)i 

The mechanism and the symmetry of pairing have 
been discussed in several theoretical studies^SiLSiS. POS 
should be distinguished from the other unconventional 
superconductors, in that it has a non-magnetic ground 
state of the localized /-electrons in the crystalline elec- 
tric field. The origin of heavy Fermion behavior in this 
compound has been discussed in terms of the interaction 
of the electric quadrupole moments of Pr 3+ with the con- 
duction electrons, rather than local magnetic moments 
as in the other heavy Fermion superconductors. There- 
fore the relation between the superconductivity and the 
orbital fluctuation of /-electron state has aroused great 
interest; POS is a candidate for the first superconduc- 
tor mediated neither by electron-phonon nor magnetic 
interactions. Hence it is of the utmost importance to 
determine the symmetry of the superconducting gap. At 
present, however, the pairing symmetry of POS is still un- 
clear. This is simply because we lack both experimental 
data and theoretical analysis enough to address the pair- 



ing symmetry. So far, a possibility of anisotropic s wave 
symmetries has been discussed for spin-singlet Cooper 
pairs. In the low-temperature phase (L-phase), an addi- 
tional symmetry breaking decreases the number of point 
nodes to 4 or 2 6,9 . The spin-triplet superconductivity 
still has a possibility, where the pairing interaction is 
mediated by the quadrupolar fluctuations. The double 
transition is more easily constructed in spin-triplet pair- 
ing with degeneracy due to the time-reversal symmetry 
than the spin-singlet pairing^. In a theorjj^, unitary and 
nonunitary spin-triplet states are proposed for H- and 
L-phase, respectively. 

Generally speaking, the tunneling spectra are expected 
to reflect the bulk density of states (DOS) of supercon- 
ductors. This is true for isotropic s wave superconduc- 
tors. In unconventional superconductors, however, the 
tunneling spectra often differ from the bulk DOS. A zero- 
bias conductance peak (ZBCP) of high-T c materials is 
an important examplo 10 ! 11 : 12 : 13 : 14 : 15 : 16 ! 17 ! 18 . The inter- 
ference effect of a quasiparticle enables the zero-energy 
Andreev bound states on the Fermi energy at surfaces 
of d wave superconductorsiSiSS. The formation of the 
zero-energy states (ZES) is a universal phenomenon ex- 
pected in unconventional superconductorsiiiSiiS 2 ^ 2 ^ 2 ^ 2 ^, 
and affects the low-temperature behavior of charge trans- 



port properties 



26,27.28,29,30 



and the Josephson current 



31,32 1 33,34 1 35,36 1 37 1 38 1 39 1 40 1 41 1 42 1 43 1 44 - When the dire ction of 

the electric current deviates from the a axis of high-T c 
superconductors, a large conductance peak is observed 
around the zero-bias, which reflects the DOS of such sur- 
face states. When the current is parallel to the a axis, on 
the other hand, the conductance shape is close to that 
of the bulk DOS in high-T c superconductors. Thus the 
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tunneling spectra are essentially anisotropic in unconven- 
tional superconductors, which means that it is possible to 
extract useful information of the pairing symmetry from 
tunneling spectra. 

The purpose of this paper is to demonstrate the dif- 
ferential conductance in normal-metal / insulator / POS 
junctions for several pair potentials proposed in recent 
studies. The junctions are described by the Bogoliubov- 
de Gennes equation^ 5 , and the conductance is calculated 
from the normal and the Andree^^ reflection coefficients 
of junctions. We discuss candidates of pair potentials in 
anisotropic s wave symmetry for the spin-singlet pairing. 
The conductance is sensitive to the relation between the 
directions of currents and the position of point nodes. 
In some cases, shapes of the conductance deviate from 
those of the bulk DOS. In the spin-singlet pairing, we 
found that the conductance vanishes in the limit of the 
zero-bias for most candidates. While in the spin-triplet 
pairing, we discuss the conductance for several candi- 
dates of pair potentials in H-phase and in L-phasc. The 
results show peak structures in the sub gap conductance 
for all candidates. 



II. MODEL 

We consider a junction between a normal metal (left 
hand side) and a POS (right hand side) as shown in 
Fig. ^ The geometry is chosen so that the current flows 
in the z-direction. Periodic boundary conditions are as- 
sumed in the transverse directions to the current and the 
cross section of the junction is S. The junction is de- 
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FIG. 1: The normal metal / POS junction is schematically 
illustrated. 



scribed by the Bogoliubov-de Gennes (BdG) equation 45 , 
fdr' 

= E 
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u{r) 
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In POS, the pair potential is expressed in the Fourier 
transformation 



A(H,r- r )=i^A(fe)e 



A(k) = 



id(k) ■ &&2 '■ triplet 



id(k)&2 



singlet, 



(3) 
(4) 



for Z c > L, where L is the thickness of the insulator 
as shown in FigH R = (X C ,Y C ,Z C ) = (r + r')/2, and 
r r = r — r' . The anisotropy of the pairing symmetry is 
characterized by A(fc) with k = (p, k z ) and p = (k x , k y ). 
In normal metals and insulators, the pair potential is 
taken to be zero. The unit matrix and the Pauli matrices 
are denoted as <To and &j with j = 1,2,3, respectively. 
Throughout this paper, we measure the energy and the 
length in units of the Fermi energy, fj,p = h 2 kp/2m, and 
respectively. The potential of the insulator is given 

by 



V(r) = V [Q(z) - Q(z - L)] 



(5) 



and q z = kp^^s/ [If) — (k z /kp) 2 is the wave number 
in the z direction at the insulator. The Andreev and the 
normal reflection coefficients of the junction are calcu- 
lated analytically 



Tee = - ZqZi 



The = 
#(±) = 

A,,± = 
K,± = 



(T - W 



e-^Ak 2 z q 2 z k\ +) R (+) | 21 |^ - z 2 W 



\ Zl \ 2 a -z 2 W 



Pi 



i..± 



i±\ 2 - (-l) l \q±l 



E*-A?.-E, 



Pl,± =\q±\fia - (-1) q± ■ <t, 

q ± ~id± x d* ± , 
z sinh(q z L), 

z\ —{q 2 - k z ) sinh(q z L) + 2ik z q z cosh(g z L), 



(6) 
, (7) 

(8) 

(9) 
(10) 

(11) 
(12) 

(13) 
(14) 
(15) 
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where ip s is a macroscopic phase of superconductor, 
k z = k z /kp, q z — q z /kp, and /(= 1 or 2) indicates 
the spin branch of a Cooper pair. These coefficients are 
characterized by the two Fourier components of the pair 
potentials, 



A 



(±) 



id±&2 : d±=d(p,±k z ) : singlet 
id± ■ <r<7 2 : d± = d(p, ±k z ) : triplet. 



In unitary states, we find 



\D±\ = 



P±\ 



E 



-co, 



\d±\ 
\d±\ 



singlet 
triplet. 



(16) 

(17) 
(18) 



The differential conductance is given b; 

2 p2tt i-n/2 

G NS (E) = —N c d<p d6 sinO 
h Jo Jo 



,47.48 



x Tr 



co 



f ee rt e + r he fl e 



B=eV b , 



(19) 



where k x = kp sin 9 cos 0, k y — kp sin 9 sin <f>, k z = 
kp cos 9, N c — Skp/(2ir) is the number of the propa- 
gating channels on the Fermi surface and Vuas is the 
bias voltage applied to the junctions. The normal con- 
ductance of the junction is also calculated to be 



2e 2 

Gn =—t-N c Tb, 
h 



T B = 



2tt 



tt/2 

d9 sin 9- 



(20) 
(21) 



, jo 4fc*# + zg' 

where Tb is the transmission probability of the junctions. 

III. SPIN-SINGLET 

Several candidates of pair potential are proposed the- 
oretically for the spin-singlet superconductivity^. Here 
we show two sets of pair potentials discussed in Ref. 0, 



d(Hl) =A - (1 - hi - k* 
d(Ll) =A (1 - k* - ki) 



d(H2) 
d{L2) 



=A (1 
=A (1 



(22) 

(23) 
(24) 
(25) 



where Ao is the amplitude of the pair potential at the 
zero temperature, kj = kj/kp for j = x,y and z are 
the normalized wave numbers on the isotropic Fermi sur- 
face. When H-phase is described by d(Hl) (d(H2)), cor- 
responding L-phase is characterized by d(Ll) (d(L2)). 
In these pair potentials, anisotropic s wave symmetry 
is assumed to have a number of point nodes. In what 
follows, we define 'node directions (n,^)' in which the 
pair potential has point nodes. The pair potential of 



d(Hl), for instance, has 6 point nodes. The node direc- 
tions are (k x , k y , k z ) = (±1, 0, 0), (0,±1,0), and (0,0, ±1). 
The thermal conductivity experiment indicates at least 6 
point nodes in H-phase. In Fig. El we show the tunneling 
conductance of the pair potentials for d{Hl) and d(Ll) 
for several choices of Lkp. Throughout this paper, we 
fix Vo/hf = 2.0 and choose three values of Lkp such as 
0.0, 0.5 and 2.0. The transmission probability of junc- 
tions, T s , are about 1.0, 0.4 and 0.003 for Lkp = 0.0, 0.5 
and 2.0, respectively. In what follows, the junction with 
Lkp = 2.0 is refereed to as the low transparent junction 
or the junction with Tb <C 1. The results in Fig. [21 (a) 
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FIG. 2: The tunneling spectra of d(Hl) in (a) and those of 
d(Ll) in (b)-(c). In (b), the current is parallel to the node 
directions of d(Ll). In (c), the current is perpendicular to the 
node directions of d(Ll). The transmission probability of the 
junction in the normal states is denoted by Tb- 



are the conductance for Eq. (1221 . In the limit of Tb <C 1 , 
the conductance shape is close to that of the bulk DOS 
denoted by a dot-dash line. Here the density of states 
are normalized by those of the normal state at the Fermi 
energy. When the pair potential are given in Eq. 12311 . 
the conductance depends on the current direction. In 
Fig. [3(b), the current is parallel to the node direction of 
d(Ll) (i.e., I j /n n d). The conductance shape in the limit 
of Tb <C 1 becomes similar to that found in Fig. [21 (a). 
When the current is perpendicular to the node direction 
(i.e., I-Lrind), on the other hand, the large enhancement 
of the conductance is seen at E = Ao as shown in Fig. [21 
(c). Thus the tunneling spectra become anisotropic be- 
cause of the anisotropy in the pair potential. The con- 
ductance shapes deviate from those of the bulk DOS even 
in the limit of T B < 1 in both Figs. [21 (b) and (c). 

In Fig. [31 we show the tunneling spectra for Eqs. I|24|) 
and l|25|l . The pair potential of d(H2) is equivalent to 
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d(Ll) under an appropriate rotation. Thus Fig. |3] (a) 
and (b) are the same with Fig.|2(b) and (c), respectively. 
There are two point nodes in the direction of (0,±1,0) in 
d(L2). In Figs. 13(c) and (d), the current is parallel and 
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FIG. 4: The tunneling spectra of d(L3) are shown in 8a) and 

(b) . The current is perpendicular to the node directions in (a). 
In (b), the current is perpendicular to the node directions. In 

(c) , the conductance is plotted for d(H3). 



FIG. 3: The tunneling spectra of d(H2) in (a)-(b) and those 
of d(L2) in (c)-(d). In (a) and (c), the current is parallel to the 
node directions. In (b) and (d), the current is perpendicular 
to the node directions. 

perpendicular to the node directions of d(L2), respec- 
tively. In Fig. |31 (d) , there is a large peak at E = Ao and 
the sub gap conductance has the U-shape as that of the 
bulk DOS. On the other hand in (c), the singularity at 
E = A is slightly suppressed and the sub gap conduc- 
tance has V-shape. In Fig. [31 (c) and (d), the anisotropy 
of the pair potential mainly appears in the shape of the 
sub gap conductance. 

When H-phase is characterized by Eq. (|22(1 . an 
anisotropic s + id wave pair potential in L-pahse is pro- 
posed by Goryo^, 



d(L3) =A 



(26) 



In the second term, the d-wave component multiplied by 
i breaks the time-reversal symmetry. The pair potential 
in Eq. H26J) has two point nodes on the Fermi surface in 
(0,±1,0) directions. In Fig. H (a), we show the conduc- 
tance for Eq. (|26|1 . where the current is perpendicular to 
the node direction. When the current is parallel to the 
node direction, the conductance is plotted in (b), where 
k\ — k 2 x in Eq. (|2*rj|) is replaced by k 2 x — k 2 y . In both (a) and 
(b), the conductance in low transparent junctions has a 
peak around E ~ A . The conductance in (a) is almost 
zero for E < 0.75Ao and are close to the bulk DOS for 
E > 0.75Aq. On the other hand in (b), the conductance 



deviates from the bulk DOS even in the limit of Tb <C 1 
and has the V-shape subgap structure. The anisotropy 
of the pair potential appears in the shape of the sub gap 
conductance as well as those in Figs.|21(c) and (d). 

In addition to Eq. (|22|1 . it is possible to consider a 
pair potential with 6 point nodes by using gap functions 
of the cubic symmetry (Oh)~- For example, a simple 
linear combination of three d wave gap functions, 

d(H3) = A (k x k y + k y k z + k z k x ) , (27) 

has 6 point nodes. We show the conductance for d(H3) 
in Fig. 01(c). The pair potential d(H3) changes its sign on 
the Fermi surface, which is the most important difference 
between Eq. (J57JI and Eqs. - |2H1> - As a consequence, 
the conductance has the ZBCP as shown in Fig. 01 (c) 
because a relation d + ~ — gL is approximately satisfied 
for|L|~l. 



IV. SPIN-TRIPLET 

As well as the spin-singlet superconductivity, a pos- 
sibility of the spin-triplet superconductivity is also dis- 
cussed in POS 7 . Ichioka et. al. proposed a pair potential 
for H-phase^, 

rn - - - - - 

d(Hl) =A W — (k x + ik y )(ky + ik z )(k z + ik x )e u (28) 
V o 

where &\, e 2 and e 3 are three unit vectors in the spin 
space. Although Eq. (|28|1 is not included in the gap func- 
tions of cubic symmetry (Oh), it explains 6 point nodes 
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on the k x , k y and k z axes. In Fig. |SJ we show the con- 
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FIG. 5: Fhe tunneling spectra of d(Hl). The transmission 
probability of the junction in the normal states is denoted by 
T B . 

ductance for the spin-triplet pair potentials in Eq. i|28|) 
for several choices of Tg ■ When the d vector has a single 
component, Eq. (|13(1 becomes 



W 



l+\\d- 



d± =e\d ± \e 1 ' 



(29) 
(30) 



where e is a unit vector which points the direction of the d 
vector. In the case of e l *'- _t * > + = — 1, the ZBCP appears 
because of the ZES^. When e^-"^+ = 1, on the other 
hand, a peak- like structure is expected around E = Aq. 
In Eq. (|28|1 . e l ^~ -l ^ + is a complex value because the pair 
potential breaks the time-reversal symmetry. In such a 
situation, the resonance energy deviates from both E = 
and E = Aq and the resonance peak is expected between 
E = and E = A(jSi. As a result, the conductance peak 
can be seen in the sub gap region as shown in Fig. [5] The 
bulk DOS vanishes at E = 0, whereas the conductance 
remains a hnite value even in the limit of Tg <C 1 , which 
reflects the surface states due to the interference effect of 
a quasiparticle. 

When H-phase is described by Eq. lj2l31) , corresponding- 
pair potential in the L-phase are given by 



d(Ll) =A [{k x + iky)(k y + ik z )(k z + ik x )ei 

+~k x e 2 ] , (31) 



d(L2) =A [(k x + iky)(k y + ik z )(k z + ik x )e x 

+(k x + ih)e 2 ]. (32) 



FIG. 6: The tunneling spectra of d(Ll) in (a)-(b) and those 
of d(L2) in (c)-(d). The current flows in the node direction in 
(a) and (d). In (b) and (c), the node direction is perpendicular 
to the current. 



These pair potentials are in the nonunitary states. In 
L-phase, some of point nodes are removed by adding the 
p wave component to the d vector in Eq. I|28|) . There are 
4 and 2 point nodes in Eas. l|31|l and (|32|l . respectively. 
Since it is difficult to determine the relative amplitudes 
of ei and e 2 components, we simply add them with an 
equal amplitude. In Fig. [SI we show the conductance in 
these L-phase pair potentials. When the current flows 
in the node direction of Eq. I]3ip. the results are plot- 
ted in (a). The conductance for small Tg has a peak 
around E = 0.3Ao which may come from the large peak 
in Fig.[SJ The DOS has a small peak at E — 1.3Ao which 
corresponds to the maximum value of Ai ,± in Eq. @- In 
(b), the current is perpendicular to the node direction of 
Eq. (|3 1 p> . where e 2 component in Ea. H31|) is replaced by 
k z e 2 . The conductance for small Tg has a large ampli- 
tude around the zero-bias. In spin-triplet superconduc- 
tors, d = — d + represents the condition for the perfect 
formation of the ZES. Actually when d + = d = vd- 
with v = ±1, the Andreev reflection probability becomes 



R A = Trr he r j he = ^ 



l 



4Pg 2 A z X, 



4Pg 2 A2+ z 2( A 2_^ 



(33) 



In the limit of E — > and zq 1, this goes to 

Ra = { V ° > (34) 



-1, 
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where spin degree of freedom give rise to a factor 2. Thus 
the zero-bias conductance is independent of Tg when 
d = d + is satisfied. The pair potential in (b) par- 
tially satisfies the condition because the e 2 component is 
an odd function of k z . As a consequence, the conductance 
at E = increases with decreasing Tb as shown in (b). 
Thus the anisotropy of the pair potential in Eq. I|31[) ap- 
pears the conductance shape around the zero-bias. The 
conductance for Eq. ij3"^|) has a large peak as shown in 
(c), which is also explained by the ZES. On the other 
hand, the conductance linearly decreases with decreas- 
ing E in (d), where ik x + ik z )e 2 in Eq. (|32[1 is replaced 
by (k y + ik x )e 3 . A peak around E = O.8A0 may come 
from the large sub gap peak in Fig. We note that the 
position of the sub gap peaks may depends on parame- 
ters such as the thickness of the insulating layer and the 
relative amplitudes among the components in d vectors. 
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FIG. 7: The tunneling spectra of d(H2) in (a), d(H3) in (c) 
and d(H4) in (d). The conductance for BW states is shown 
in (b) for comparison. 

For H-phase, there are another candidates of the pair 
potentials such as^ 

d(H2) =A [k x ei + k y ee 2 + k z e 2 e 3 ] , (35) 
d(H3) =2A [k x {P z - P y )e x + k y {P x - P z )e 2 

+ k z (P y -P x )e 3 ] ) (36) 
d(HA) =2A [k x (k 2 - k 2 y)ei + k y {P x - P z )ee 2 

+ k z Ck 2 y -k 2 x )e 2 e 3 }, (37) 

where e = e l27r / 3 . The pair potential in Eq. (|35|l is similar 
to that of Barian-Werthamer (BW) states^ described by 

d(BW) = A [k x ei + k y e 2 + k z e 3 ] . (38) 



Eq. (|35|l . however, is in the nonunitary states because of 
a phase factor. One spin branch has a full gap, other 
has 8 point nodes in (±1,±1,±1) directions. The node 
directions of this pair potential contradict to the exper- 
imental results. In Fig. [7| we show the conductance for 
Eqs. H35JI in (a). For comparison, we also show the con- 
ductance of the BW states in (b). The conductance for 
Tb -C 1 increases rapidly with increasing E and has 
a peak around E = 0.2Ao as shown in (a). We note 
that the conductance spectra of the BW state in (b) also 
show a peak around E = 0.5Ao- The peak structure 
may indicate some surface states of the BW type super- 
conductors because the bulk DOS only have a peak at 
E = max(Ai,±) = 1.4A in (a) and E = A in (b). 
When d vectors have more than two components, the 
shapes of the conductance spectra tend to have sub gap 
peaks. Mathematically speaking, when d_ is not parallel 
to d* + , the product of the two pair potentials in Eq. H13|) 
becomes 

A ( _)At +) = d ■ d* + a + id- x d* + ■ <x. (39) 

The second term is a source of the sub gap peaks in the 
BW type states. At present, however, we have not yet 
confirmed an existence of some surface states. In Fig. [7| 
(c) and (d), we show the conductance for Eqs. il3T)|) and 
(|37|l . respectively. There are 14 point nodes on the Fermi 
surface in Eqs. (|36() and (|37|l . Although the number of 
point nodes are larger than that found in the experi- 
ment, these pair potentials explain the 6 point nodes in 
k x , k y and k z directions. The conductance in Fig. 0(c) 
shows peak structures at E = 0.36Ao and 0.76Ao- These 
peaks are far from a peak in the bulk DOS at E = Aq. 
The conductance in Fig. [7| (d) also shows peak structures 
at E = 0.13A , 0.34A and 0.86A . However, there is 
no structure in the bulk DOS around the lowest peak. 
In addition to Eqs. (|35|I - H37|) . the polar state and the 
Anderson-Brinkman-Morel (ABM)£i state are proposed 
for H-phase of the spin-triplet pairing 7 , 

d(polar) =A k z e 3 , (40) 
d(ABM) =A (k x + ik y )e 3 . (41) 

The transition to L-phase is caused by the spin-orbit cou- 
plingi. The polar state in Eq. (|4U|I has a line node at 
k z = and the ABM state in Eq. (|41|> has two point 
nodes at k z = 1. In Fig. [S] (a), we show the conduc- 
tance in Eq. (|40|l . where a plain including the line node, 
k z = 0, is perpendicular to the current. The results show 
the ZBCP because Eq. l|4T?|) satisfies d = d + . In (b), 
we show the conductance in the ploar state, where k z in 
Eq. d40|) is replaced by k x and a plain including the line 
node, k x = 0, is parallel to the current. The conduc- 
tance at the zero-bias vanishes in the limit of Tb <C 1 
and increases linearly with increasing E. The shape of 
the conductance, however, deviates from that of the bulk 
DOS. In Fig. |S] (c), we show the conductance in Eq. (|4*T1 . 
where the node direction is parallel to the current. In 
low transparent junctions, the conductance vanishes in 
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FIG. 8: The tunneling spectra of d(polar) in (a) and (b). 
Those for d(ABM) are in (c) and (d). 

the limit of E — ► 0. The shape of the conductance, how- 
ever, deviates from that of the bulk DOS. In (d), we 
show the conductance in the ABM state, where k x + ik y 
in Eq.l^TJ is replaced by k z + ik x and the node direc- 
tion is perpendicular to the current. The broad ZBCP 
appears because Eq. (|41(l satisfies d_ = — d+ only when 
\k z \ = l 23 . The height of the ZBCP is expected to be 
much larger in junctions with thicker insulating layers. 
The transmission probability for perpendicular injection 
to the thicker insulating layers become much larger than 
those for another incident angles. As a consequence, the 
condition d = d + is better satisfied in junctions with 
thicker insulators. 

In comparison with the spin-singlet pairing, the con- 
ductance in the spin-triplet superconductivity tends to 
have the sub gap structures. The peak structures in 
Figs. [SJ IH1 (a) and El (d) are stemming from the broken 
time-reversal symmetry states in Eq. <|28|) . The ZES is 
responsible for the peaks around the zero-bias in Figs. El 
(b) and (c). The d vectors with multi components are the 
origin of the peaks in Fig. El Thus POS may be the spin- 
triplet superconductors if the sub gap conductance shows 
complicated peak structures in experiments. The argu- 
ment, however, is still a guess based on the calculated 
results. This is because it may be possible to consider 
another pair potentials with 6 point nodes. 

In this paper, we do not consider the self-consistency 
of the pair potential near the junction interface. It is 



empirically known that the depletion of the pair potential 
modifies the conductance structure around E = Ao or 
maximum of Ai ±. Our conclusions remain unchanged 
even in the self-consistent pair potential unless the self- 
consistency does not change the symmetry of the pair 
potential and/or the number of components in d vectors. 



V. CONCLUSION 

We have discussed the differential conductance in 
normal-metal / insulator / POS junctions based on the 
Bogoliubov-de Gennes equation. For spin-singlet pairing, 
the conductance is calculated for three candidates of pair 
potentials in the anisotropic s wave symmetry. The re- 
sults show that the conductance spectra depend strongly 
on the relation between the direction of currents and that 
of nodes. We found that the conductance vanishes in the 
limit of the zero-bias and there is no anomalous behavior 
around the zero-bias for these candidates. The conduc- 
tance for s + id wave symmetry in L-phasc and that for 
d wave symmetry in H-phase are also demonstrated. In 
the case of spin-triplet superconductivity, we discuss the 
conductance for six candidates of pair potentials in bi- 
phase and two candidates in L-phase. The results show 
peak structures in the sub gap conductance for all can- 
didates. The broken-time reversal symmetry states, the 
zero-energy states and d vectors with multi components 
arise these peak structures. POS may be a spin-triplet 
superconductor if the peak structures in the sub gap con- 
ductance is observed in future experiments. In particular, 
the presence or the absence of the ZBCP is an important 
information to address the pairing symmetry. Thus ex- 
periments of the tunneling spectra are desired^. 

Recently, it has been pointed out that the magneto 
tunneling spectroscopy is a useful tool to know details of 
internal structures of pair potentials 5 -. The tunnel spec- 
tra through a ferromagnetic tip 5 * reflect the spin con- 
figuration of Cooper pairs in the case of the spin-triplet 
superconductors. Even in the spin-singlet superconduc- 
tors, the absence of the time-reversal symmetry in ferro- 
magnets affects the interference effects of a quasiparticle 
and modifies the tunneling spectra. At present, however, 
the investigations in this direction are limited in the d 
wave high-T c superconductors. The magneto tunneling 
spectroscopy in another unconventional superconductors 
is a future problem. 

In this paper, we assumed the clean ballistic junctions. 
It is known impurity scatterings in normal metals induce 
the proximity effect^ 5 -. The proximity effect of unconven- 
tional superconductors with point nodes is also a open 
question. 
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